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The Perturbative Approach (PA) introduced by lAlardl (120071) provides analytic solutions 
for gravitational arcs by solving the lens equation linearized around the Einstein ring solu- 
tion. This is a powerful method for lens inversion and simulations in that it can be used, in 
principle, for generic lens models. In this paper we aim to quantify the domain of validity of 
this method for three quantities derived from the linearized mapping: caustics, critical curves, 
' and the deformation cross section (i.e. the arc cross section in the infinitesimal circular source 

approximation). We consider lens models with elliptical potentials, in particular the Singular 
Isothermal Elliptic Potential and Pseudo-Elliptical Navarro-Frenk- White models. We show 
that the PA is exact for this first model. For the second, we obtain constraints on the model 
parameter space (given by the potential ellipticity parameter e and characteristic convergence 
k s ) such that the PA is accurate for the aforementioned quantities. In this process we obtain 
analytic expressions for several lensing functions, which are valid for the PA in general. The 
IT^i ■ determination of this domain of validity could have significant implications for the use of the 

PA, but it still needs to be probed with extended sources. 
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1 INTRODUCTION 

Gravitational arc systems can be used as unique probes of the 
matter distribution of galaxies and galaxy clusters acting a s 
lenses jKovnei|[l989l : iMiralda-EscudeJI 1993ar . lHattori et all 1997b . 
Further, the ir abundance can be used to constrain cosmologi- 
cal models dBartelmann et alj|l998l : lOguri et afooil ; iGolse et all 
|2002| ; iBartelmann et al. 2003). This motivated several arc searches 



to be carried out, both i n wide field survey s (IGladders et al 
20031: lEstrada et al.l l2007t ICabanac et al.l |2007|; iBelokurov et al 



20091: iKubo et alj l201tj: iKneib et al] l20icl: iGilbank et alj 1201 ll; 
Wen et al.l |201 ll: More et al.l 120121: iBavlissI [20121; IWiesner et al] 



2012l:lErben et al.ll2012h. as well as in images targeting know clus- 
ters dLuppino et aLlll999b IZaritskv & Gonzalezl 120031; ISmith et all 



20051: ISand et alj 120051: iHennawi et al]|2008l : iKausch et alj|20ld : 
Horesh et all 120101 :1 Furlanetto et al.l 120120 . Upcoming wide field 



imaging surveys , such as the Da rk Energy SurvejQ (DES; 
lAnnis et al]|2005l : lAbbott et al]|2005h . which started operations in 
2012, are expected to detect of the order of 10 3 strong lensing sys- 
tems, about an order of magnitude increase with respect to the cur- 
rent largest surveys. 

Two primary approaches have been followed in practical 
applications of gravitational arc systems. On the one hand, in- 
verse modelling attempts to "deproject" the arcs i n individual 



19931; iKeetonl l200ld: IGolse et alj |2002|; 


Comerford et alj 120061: 


Wavth & Webstei l2006t Julio et alj 2007 


, l2010h. On the other 


hand, arc statistics (Wu & Hammer 1993; Grossman & Saha 1994; 
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iBartelmann & W eiss 1994) aims at counting the number of arcs in 
cluster samples and comparing with the predictions from cosmo- 
logical models. 

Both approaches require the lens equation to be solved nu- 
merically for finite sources numerous times. The inverse modelling 
typically needs arc images obtained from a multidimensional space 
of source positions and lens parameters scanned during the min- 
imiza tion process to find the best solution for sources and lenses 
(e.g., IWavth & Webste3l2006h . For this reason analyses using the 
inverse modelling are often restricted to simple lens models, in 
particular models with elliptic lens potentials (so-called pseudo- 
elliptical) and/or to point sources, for example, co nsidering bright 
spots in arcs as multiple images of point sourc es jKeetonll200M 
IWavth & Websterl2006l : ljullo et al.ll2007l ; lOgurfcOlOl) . 

For arc statistics, predictions for the arc cross section 
must be obtained as a function of source and lens properties 



arc images dMiralda-Escude 


Il993bl; IBartelmann & Weissl Il994l; 


iMeneghetti et al. 200 ll 12003 


: Oguri et al .1 2003). The cross sec- 



tion is then convolved with the distribution of lens properties ex- 
pected in a given cosmology and convolved with the source dis- 
tribution. Another approach is to use directly high resolution N- 
body simulations obtaining arc images by ray- tracing through the 
mass distribution for a l arge number of sour ces dHoresh et al .120 111 ; 
IMeneghetti et"aT]|2008l : iBoldrin et al]|2012h . 

It is therefore useful to develop approximate methods for ob- 
taining gravitational arcs, which will be particularly useful given 
the increase of strong lensing systems to be discovered by the next 
generation wide-field surveys. A most promisin g technique for thi s 
purpose is given by the Perturbative Approach ( Alard l2007ll2008l) . 
which provides an approximate solution for the lens equation close 
to the Einstein ring, leading to analytic solutions for arcs. The 
power of this approach is that it can be applied, in principle, to 
generic lens models, including those arising from simulations. The 
method is suitable for large tangential arcs, since the solutions are 
accurate for images located close to the Einstein ring corresponding 
to the circularly averaged lensing potential. 

Another important feature of the method is that it naturally re- 
produces arcs resulting from the merger of multiple images, which 
cannot be accounted for with o ther approximate methods for arc s 
proposed in the literature (e.g., iKeetorfeOOlbl : iFedeli et al.ll2006h . 
Such merger arcs are key for lens inversi on methods and a lso play 
an important role in the arc cross section dRozo et al.ll2008l) . 

The Per turbative Approach has already been used for inverse 
modelling in[Xlard (2009, 201(1). Given that it reproduces arc con- 
tours that can be associated to isophotes, it could also be used to 
simulate the brightness di stribution of arc s , in a similar way to 
what was implemented in iFurlanetto et all I d2013l) for arc shaped 
contours. 

An important issue for practical applications of this approach 
is the det e rmina tion of its domain of validity. This topic is discussed 
in lAlardl d2007h . com parisons with arc sim ulations are presented 
inlPeirani et al.l J2008b . and a recent work bv lHabara & Yamamotol 
d201ll) has investigated arcs in several configurations for a pseudo- 
elliptical model in this approach. However a systematic study of 
its limit of applicability has not yet been carried out. In this paper 
we make a first attempt to determine a domain of validity of the 
method in terms of the parameter space of the lens model. We will 
restrict to the simple case of pseudo-elliptical models, which are 
nevertheless widely used for the inverse modelling. Moreover, for 
simplicity, we will restrict the comparisons with the exact solution 
for three quantities connected to arcs, but which do not involve the 



lensing of finite sources. We expect that the limits here obtained 
can be connected to the domain of validity for arcs and extended 
for more general models, but this is left to subsequent explorations. 

In this work, our purpose is twofold. The first is to explore the 
application of the Perturbative Approach to determine quantities 
arising from the local lens mapping, such as the arc cross section 
for infinitesimal circular sources (deformation cross section). The 
second is to determine a domain of validity such that the critical 
curves, caustics, and deformation cross section are accurately ob- 
tained. This study is performed for the pseudo-elliptical Navarro- 
Frenk-White model (PNFW), determining regions of its parameter 
space where the Perturbative Approach provides a good approxi- 
mation for these quantities. We also consider the Singular Isother- 
mal Elliptic Potential (SIEP) model and show that the solution of 
the Perturbative Approach is exact in this case. 

The outline of this paper is as follows: in Sec. [2] we present a 
few basic results of gravitational lensing theory, introduce the radial 
lens models to be used in this work, and discuss models with elliptic 
lensing potentials. In Sec. [3] we review the Perturbative Approach, 
present its application to the computation of the deformation cross 
section, and discuss its implementation in pseudo-elliptical mod- 
els. In Sec.|4]we establish a metric for the comparison between the 
Perturbative Approach and the exact solution for critical curves and 
caustics and determine a domain of validity for the Perturbative ap- 
proach. In Sec. [5] we summarize the results and present concluding 
remarks. 



2 BASICS OF GRAVITATIONAL LENSING: 
DEFINITIONS AND NOTATION 

In this section we present a brief review of the lensing theory to 
set up the notation and to define the quantities associated with 
pseudo-elliptica l mode ls . For a more detailed descri pt ion see, e.g., 
ISchneider et al. | dl992l) : iMollerach & Rouletl d2002l) ; IPetters et al.l 
d200ll) . 



2.1 Lens Equations 

The standard description of strong lensing phenomena is based on 
the lens equation, which relates the transverse position (with re- 
spect to the optical axis) of the observed images £ to those of the 
sources rj, where the boldface denotes two-dimensional vectors, in 
the image and source plane, respectively. We may choose a length 
scale £o and define x = £/£o and y — 7]/r]o, with 

Dos. 



'/<> 



£>OL 



6>, 



where Dol , -Dos are the angular diameter distances from the ob- 
server to the lens and source respectively. Using these definitions 
the lens equation is written as 



y = x - a(x) = x - V^ipix), 



(1) 



where a(x) is the "dimensionless" deflection angle and <p(x) is 
the "dimensionless" lensing potential, 

The local distortion in the lens plane is described by the Jaco- 
bian matrix of eq. |[T) 



dy 
Ox 



8 t j - diOtj(x) 



(2) 



The two eigenvalues of this matrix are written as A r = 1 — k+7 and 
At = 1 — k — 7, where k and 7 are the convergence and the shear 



On the validity of the Perturbative Approach for Strong Lensing 3 



given below. Magnifications in the radial and tangential directions 
are given by the inverses of these eigenvalues. Points satisfying the 
conditions X r ,t = define the radial and tangential critical curves 
respectively. Mapping these curves onto the source plane, we obtain 
the caustics. 

For axially symmetric models the lensing potential is a radial 
function, i.e., y>o(x) = ipo(x). Hence, the deflection angle is given 
by 



t(aj) = 



dx 



(3) 



where ~£(£ox) is the mean surface density within a radius x and 
E cr it is the critical surface density defined by 



Ecrit 



.Dos 



4ttG D ls D ol 

with Dls being the angular diameter distance between the lens and 
the source. The convergence and shear can be written in terms of 
the deflection angle as 



a(x) da(x) 



dx 



y(x) 



a(x) da(x) 



dx 



(4) 



The tangential critical curve (i.e., the "Einstein Ring" in this case) 
is a circle with radius xe, obtained by solving 



Xt(x) 



y(x) 



0. 



(5) 



which is a function of the parameter k s alone. The deflection angle 
and eigenvalues \,. : t can be straightforwardly derived using eqs. ((3} 



2.2 Pseudo-Elliptical Models 

Models with elliptical potentials p rovide simple analytica l solu- 
tions for some lensing quantities jBl andford & Kochanefc! 1 19871 : 
iKassiola & KovneJ 1 19931 ; iKneitJ l200ll) . They have been widely 
used in lens inversion problems and are implem ented in severa l 
public cod es for lens inversion such as Gravle ns dKeetonl 200 It ) 



Lensview dWavth & Webster! 120061) , Lenstool l ljullo et al 
and glafie (Oguri 2010). They have also been used for 



2007) 



arc simu- 



lations dOgur i 2002; Me neghetti et alj200ll2007l) . 

The standard prescription for obtaining an elliptical potential 
<p E (x) from a given axially symmetric one tpo (x) is to replace the 
radial coordinate x in the latter by 



ai x\ + a,2 x\ = 



where 

= \/ 

such that the ellipticity of the lensing potential is 



(11) 



(12) 



The radial critical curve is determined by 

da(x) 



1 



dx 



0. 



(6) 



Going to the source plane, the tangential caustic degenerates into a 
single point, while the radial caustic becomes a circle. 

In this work, one model we will make use of is the Singu- 
lar Isothermal Sphere (SIS), which is useful to model lenses at 
the galactic scale. Its dimensionless lensing potential is given b y 
dTurner et all 1984 ISchneider et~ai]l992t Ivan de Ven et al"l2009h 



<fio(x) 



(7) 



where we choose the Einstein Radius to be the characteristic scale 

2 



Co = Re = 



GEci-it 



where a v is the velocity dispersion. The deflection angle, conver- 
gence, and shear for this case are 



a(x) = 1, 



k(x) = — 
v ; 2x 



7 (a;). 



(8) 



From these expressions the eigenvalues A r , t are easily derived, 
as well as the solutions of the lens equations for finite sources 
dDobler & Keetoijl200d :l Montovj[20Tll) . 

We will also make us e of the Navarro-Frenk-White model 
dNavarro et ail 1 19961 1 19971 hereafter NFW), often used to repre- 
sent lenses in the galaxy to galaxy cluster mass scales. This model 
has two independent parameters r s and p s . By fixing £o = r 3 and 
defining the characteristic convergence as 



E c rii 



the lensing potential is given by dBartelmannl 19961) 



<fo(x) — 4k s 



1, 2 X 

2 l0g 2 



2 arctanh^ 



1 + x 



(9) 



(10) 



where the orientation was chosen such that the major axis of the 
ellipse is along the x 1 axis (i.e., a,2 > ai). 

The deflection angle, convergence, and shear can be writ- 
ten as combinations of the lensing functions of the con'espond- 
ing axially symmetric mo del for any choice of ai and a2 
dDumet-Montova et alj2012T ). For example, the SIEP model, which 
is obtained by introducing ellipticity in the SIS potential eq. (0, has 
the dimensionless deflection angle given by 



a e (x) — — — (ai cos 4>, 02 sin cj>) . 

A<4 



(13) 



The PNFW model is similarly obtained by introducing ellip- 
ticity into the NFW lens potential eq. # 10j . 



3 PERTURBATIVE APPROACH 

For a given lens model, the Perturbative Approach allows one to 
obtain analytic solutions for arcs as perturbations of the Einstein 
Ring solution. In this work we investigate the limits of applicability 
of the Perturbative Approach, by considering simple non-axially 
symmetric models and by looking at local properties of the lens 
mapping, instead of lensed finite sources. 

In this section we briefly review the Perturbative Approach 
and use it for the derivation of the caustics and critical curves, the 
deformation cross section and quantities needed for its computa- 
tion. The method is also applied to models with elliptical lensing 
potentials. 



3.1 Lens Equation 

The gist of the Perturbative Approach for gravitational arcs devel- 
oped bv lAlardld2007ll2008h is to obtain an analytic solution for the 
lens equation considering the lens as a perturbation of an axially 
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0. 



(14) 



symmetric configuration and the source position as a small devia- 
tion from the optical axis (i.e., positioned transversely away from 
perfect observer-lens-source alignment). In other words, the arcs 
are found as perturbations of the Einstein Ring configuration. 

The Einstein Ring is the image of a source aligned with an 
axially symmetric lens (with lensing potential ipo). Its radius xe is 
obtained by solving eq. l[5} at the centre of the source plane, i.e. 

dip a 
dx 

Arcs can be obtained by perturbing the equation above either by 
shifting the position of the source away from the optical axis 
and/or by adding a non-circular perturbation to the lensing poten- 
tial. These perturbations are described by 

V = Sy, f(x) = ip (x) + 5ip(x), (15) 

and they are assumed to be of the same order throughout the fol- 
lowing calculations. The response to such perturbations is given by 
the displacement of the radial coordinate in the lens planqj, i.e., 

x = ie — > x = xe + Sx. (16) 

To find 8x we solve eq. l[T) by expanding the solution around x = 
xe- Expanding the lensing potential in a Taylor series around x = 
ie, we have 



<p(x)=J2[Cn+f n (cf>)](x-XE) n , 



where 



1 



d"<p 



dx r ' 



1 



d n Si> 



(17) 



(18) 



Inserting eqs. d!6t and d!7t into eq. (TJl, we find that the resulting 
equation at zeroth order in the perturbations is 



X'e — Ci, 



(19) 



which is the Einstein Ring equation. Using this relation we find the 
equation at the first order as 

2/i = (k 2 <5;e - /x) cos0 + — ^ sin0, (20) 
xe dip 

2/2 = (k2Sx — /i)sin</> — —7 cos(j>, 

xe dcp 

where K2 = 1 — 2C2. From eqs. ((3) and© we have 

d 2 <po(x) _ a{x) 

dx z x 
and therefore K2 can be expressed as 

k 2 =2-2«(x B ). (21) 

Eq. ( 120b is the lens equation in the Perturbative Approach. It 
can be solved for 8x for each angular position cp of the source, 
given a perturbation described by f n (4>)- To obtain the images of a 
finite source, we must first parametrize its boundary. Then, by vary- 
ing <j> from to 2-7T, each point of that boundary is mapped to the 
lens plane through eqs. d20> . As a result, a new equation with sep- 
arated radial and angular components is formed, whose solution is 



Note that in lAlardH2007ll2008l) io = xe was used as a characteristic 
scale. This choice is equivalent to setting xe = 1 in our equations. In this 
work we have made the choice of keeping xe explicitly in the equations for 
more generality, allowing us, for example, to choose another characteristic 
scale of the problem. 



obtained straightforwa rdly dAlardll2007l 120081: IPeirani et alj|2008t 
lDumet-Montoyj|201 lb 

It is important emphasize that the solutions (xi, (pi) of eq. d20| > 
are valid only to first order in the perturbations in eq. (1151 . i.e. only 
for points near the Einstein Ring. For points far from this curve, the 
solutions are not expected to be highly accurate. For this reason, 
the Perturbative Approach is particularly useful for applications in- 
volving tangential arcs. In this work, instead of using finite sources, 
we focus on the potential applicability of this method to quantities 
based on the local mapping as a first step to quantify the differences 
with the exact solutions. 



3.2 Local Mapping 

The Jacobian matrix for the lens mapping is 



J — ( 7j~ j — (Js-^L,pol) ik {JL,pol->cart, 



(22) 



where Js^L, P oi is the Jacobian of the transformation from the lens 
plane to the source plane in polar coordinates from eq. (1201 and 
*iL, P oi-icart is the standard Jacobian matrix from polar to Cartesian 
coordinates. The calculation of the eigenvalues of the lens mapping 
is then straightforward from the equation above and they are given 
by 



At = 

Therefore, the radial coordinate of the tangential critical curve is 



(23) 



x t (<f>) — x E + Sx t (<f>) = x E + — ( fi + —-77T 

K2 \ XE dip'' 

and the parametric equations of the critical curve are simply 



(24) 



Xit = x t costp, 
x 2t = x t sin ip. 

By inserting Sxt in eq. d20t . the parametric equations of the tan- 
gential caustic are found to be 



1 d 2 f 1 df . 

V u = — UT 008 ^^ U sm ' 

xe dtp- 1 xe dtp 



(25) 



2 /o 



xe dip 2 



xe dip 



3.3 Constant Distortion Curves 



For infinitesimal circular sources, the length-to-width ratio of arcs 
can be approximated b y the ratio of the eigenvalues of the lens map- 
ping Jacobian matrix 1 Wu & Hammed [19931 iBartelmann & Weissl 



ll994l:lHamana & Futamasell 19971) 



L_ 

W 



life ("01 



A r (as) 



A t (as) 



(26) 



Under this approximation, it is possible to define a region 
where gravitational arcs are expected to form by fixing a value 
for the threshold length to width ratio i? t h- Such region is limited 
by the curves R\ = ±i? t h (constant distortion curves). Although 
the condition d26t does not hold for merger arcs (IRozo et al.l20O8t 
lFerreiral20l"bT) . nor for large or elliptical sources, the curves defined 
above still provide a typical scale for the region of arc formation. 
In this work, we adopt the common choice 7? t h = 10 (unless ex- 
plicitly stated otherwise). We denote the radial coordinates of these 
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curves as x\. They are obtained by solving R\(x) = ±R t h, with 
A r and At given in the Perturbative Approach by eq. (1231 . It follows 
that 



x t x 



Hth + 1' 



Rx — +Rth, 



(27) 



Rx 



-Rth. 



The constant distortion curves in the lens plane are therefore self- 
similar to the tangential critical curve. The mapping of these curves 
to the source plane is done by substituting Sxx = xx — xe in 
eq. ( 120b . For instance, the curve Rx = +Rth has the following 
parametric equations 



V2X = J1 J, 1, E sin < 



J_ \( «th 



dzfa 

d4> 



^sin, 



dflL 



dfo 
d<f> 



(28) 



The parametric equations of the Rx = — -Rth curve are given 
by the expressions above with the substitution R t h — 1 — >• R t h + 1. 
There is no self-similarity between these curves and the tangential 
caustics. 



3.4 Deformation Cross Section 

As mentioned in the introduction, the arc cross section is a key in- 
gredient in arc statistics calculations. Is is defined as the effective 
area in the source plane such that sources within it will be mapped 
into images with L/W ^ Rth- The definition of this area must 
take into account the image multiplicity given the source position 
(i.e. multiply-imaged r egions are counted multiple times, see e.g., 
iMeneghetti et al .120031) . The computation of the arc cross section in 
general deman ds ray-tracing simulations, which are computation- 
ally expensive {Miralda-Escude1l993bl;lBartelmann & Weissl 19941 ; 
IMeneghetti et al.ll200ll. 120031 ; lOguri et alj|2003l) . An alternative is 
to use the infinitesimal circular source approximation, eq. J26b . 
which allows the computation of the arc cross section to be car- 
ried out directly from the local mapping from lens to source plane. 
In this case, tj_R th is computed in the lens plane by 



(see, 



d 1 



Fedeli etHll 



120061 : ICaminha et al] 

H = (A r At) -1 



(29) 

I2012L 
the 



iDumet-Montova et alj 20121) . where fi — (A r M) " is 
magnification and the integral is performed over the region of arc 
formation above the chosen threshold. The quantity an th is known 
as the dimensionless deformation cross section. 

In the Perturbative Approach, the magnification can be written 
from eq. ((23} as 



,.2 



x t {4>) 



(30) 



X — Xt((j>), X > Xt((p), 



where Xt{4>) is given in eq. d24b . Inserting the equation above in 
eq. d29t and integrating the radial coordinate within the lower and 
upper limits given in eq. ( 127) , it is straightforward to obtain 



2 i-Rthr + i 
a ^- K \\ Rt h?-i? 



Xt(<t>)d<f>. 



(3D 



For axially symmetric models (x t = Kb) the cross section is 



given simply by 



crB th = 2nK2X E 



(32) 



Note that aR tli oc R th for R t h 3> 1, as expect ed from the 
behaviour of the deformat ion cross section with i? t h dRozo et all 
l2008l : ICaminha et alj|2012l) . 



The express i on a bove is exact for the SIS model 
dBartelmann et alj 1 19951) . For other axially symmetric mod- 
els this expression is still an approximation, since the curves 
Rx = iRth are obtained approximatively. 



3.5 Perturbative Functions for Pseudo-Elliptical Models 

We write the elliptical potential as 

ip e (x) = (po(x) + [ipo(x) - ^o(ic)]) (33) 
such that the perturbed potential becomes 

Sip(x, 4>) = tpo(x) — <po(x). 

From the definitions ( 118) and using the identities ((3} and ((4), it 
follows that 



fi 

dfo 
d<p 

d 2 fo 
dcj) 2 



— o(xe) — oc(xe), 

Xe 

x 2 

— ^-a(i E )(a2 — ai) sin2<^, 
2-Xe 

— — q(xe)(i2 — Ol) COS 2(f) 
xe 



(34) 



7(*e) 



x 2 

— (d2 — ai) sin 2(j> 

XE 



where a and 7 are the deflection angle and shear of the corre- 
sponding axially symmetric lens. These expressions hold for any 
parametrization of th e lensing potential elliptici ty and for any 
pseudo-elliptical lens jDumet-Montova et alj2012T) . 

For small values of the lensing potential ellipticity, eqs. d34| l 
reduce to 



fi 

dfo 
dcf> 

d 2 fo 
d(f> 2 



ai 



(7 2 



k(xe)xe cos 2<j>, 



a 2 -ai 2 . 
xe sm lq>, 



(35) 



(CL2 — o,i)xe COS 2(f). 

From eq. J24t and the expressions above, we have 

0,2 — cii ( 2 — k(xe) 



Xt{4>) = xe 



1 



2 K 2 

and inserting this into eq d31t we get 

1 



cos 2cj> 



an 



o 2 -Rth + 1 
2tvx e - 



2 1 

K 2 + 



Oil 



(36) 



■(37) 



-D-L-'.v^fl 

Thus, for small ellipticities, the deviation with respect to the axially 
symmetric case is quadratic. 

Instead of using (12 and ai it is more intuitive to express the 
results in terms of the ellipticity of the potential. Several parameter- 
izations have been used to define the ellipticity in this context . From 
now on, we adopt the convention llBlandford & KochanetJ 1 19871 ; 
iGolse & Kneibll2002l ; lDumet-Montova et alj|2012h 

1 — e, 0.2 — 1 + e, 



(38) 



where e is the potential ellipticity parameter. The connection to the 
ellipticity of the mass distribution es depends on the model. For 




Figure 1. Critical curves (top panels) and caustics (bottom panels) obtained with the Perturbative Approach (dashed lines) and the exact solution (solid lines) 
for the PNFW model: k s =0.1 and e = 0.2 (left), k b = 0.5 and e = 0.32 (middle) and k s = 1.0 and e = 0.35 (right). The values of e in the middle and 
right panels were chosen by imposing T> 2 = 5 X 10 — 4 for critical curves and caustics, respectively (see Sec.|4). 



the SIEP £e = 3 e to first order in e dKassiola & Kovnej[l993h . 
For the PNFW model this re lation depends on k s and exp ressions 
for es(e, k s ) are provided in lDumet-Montova et alj J2012h . 

Fig. Q] shows the comparison for caustics and critical curves 
between the Perturbative Approach and the exact solution for the 
PNFW model for different values of « s and s. 



3.6 Singular Isothermal Elliptic Potential 

One of the simplest and most used lens models is given by the SIEP. 
For this model, by substituting ([8} in eq. l !34t . the perturbative func- 
tions are 



fi = 
dfo 



= («2 - Ol) 



sin 2(f> 
~2A7' 



(39) 



d 2 fo 



ai a^A^ 3 — A$, 



where is given in eq. dl2t . When substituted into eqs. J20l l the 
expressions above lead to 



Vi 



(l — ^-J cos cj> and y 2 — x ^1 — sin0, 



(40) 



which are the components of the lens equation of this model with- 
out any approximation (see eqs.[T|and|13|l. Hence, the solution of 
the Perturbative Approach is exact in the case of lensing by the 
SIEP model. 



The same conclusion does not hold for the PNFW model. We 
will thus investigate the domain of validity for this model in the 
next section. 



4 LIMITS OF VALIDITY OF THE PERTURBATIVE 
APPROACH FOR THE PNFW MODEL 

Previous attempts to quantify the differences between exa ct and 
pertur bative solutions were carried out in the literature. lAlardl 
( 2007) proposed a method based on the relative importance of the 
third-order term in the Taylo r series of the gravitational potential. 
lHabara & Yamamotol d201ll) performed a qualitative analysis of a 
particular arc configuration, varying some of the system parame- 
ters and establishing criteria based on the position and multiplicity 
of the images. However, they did not define a metric to compare the 
solutions nor carry on the analysis for more general configurations. 

Investigating the domain of validity of the Perturbative Ap- 
proach with finite sources would require a large parameter space to 
be probed, including the lens and source parameters and their rela- 
tive positions. On the other hand, as a starting point, we may look 
at quantities that are dependent only on the lens, such as the tangen- 
tial caustic and critical curve and the deformation cross section (the 
latter will depend also on the choice of Rth)- Besides reducing the 
parameter space — for example, for e and n s in the PNFW case — 
it is simpler to define metrics to quantify the deviation of the per- 
turbed solution from the exact one. We expect that the constraints 



On the validity of the Perturbative Approach for Strong Lensing 7 




e e 

Figure 2. Mean weighted squared radial fractional difference T> 2 as a function of the ellipticity parameter for the PNFW model for various values of the 
characteristic convergence k s . Left Panel: For Critical Curves. Right Panel: For Caustics. 



on the domain of validity determined from the quantities above can 
be connected to those arising from the images of finite sources. 
Thus, exploring the simplest case before may provide guidance to 
the determination of the domain of validity of the method finite 
sources in the future. Setting a domain of validity from the lens 
model alone may provide a rapid method to adjudicate validity of 
the Perturbative Method a priori, just from the lensing potential, 
without the need of obtaining images of the sources. 

In this section, we shall attempt to quantify the deviation of 
critical curv es and caustics using a figure -of-merit akin to the one 
proposed in IPumet-Montova et al.l d2012h . We will then compare 
the deformation cross sections and, finally, combine the results to 
obtain limits that define a region in the parameter space of PNFW 
models where the Perturbative Approach can be used to accurately 
obtain local properties of a given lens system. 

4.1 Limits for critical curves and caustics 

To quantify the deviation of the solution of the Perturbative Ap- 
proach from the exact one for critical curves and caustics we use 
a figure-of-merit defined as the mean weighted squared fractional 
radial difference between the curves, i.e0 

V 2 ^ Ejll Wj[x E s(fa) -ZPaQQ] 2 

where XEs(4>i) and xp\{fa) are the radial coordinates of the tan- 
gential curves (either critical curves or caustics) obtained from the 
exact solution and with the Perturbative Approach, respectively. 
These are computed on a discrete set of N points defined by the 
polar angle fa. Further, Wi = fa — fa-i is a weight to account for 
a possible non-uniform distribution of points in fa 

3 Expression ^411 i s for mally equal to the one proposed in 
IPumet- Montova et alj i2012l) . where it was used to compare an iso- 
contour of k to an ellipse. Here the same expression is used to compare two 
solutions for caustics or critical curves. 



Choosing a cut-off value for T> 2 , we can define a range in e 
for which the curves obtained with both the exact and perturbative 
solutions will be similar enough to each other. The cut-off value 
is then chosen by visually comparing the exact and perturbative 
solutions for the critical curves and caustics associated with several 
values of X> 2 , for combinations of the PNFW lens parameters. 

Before presenting the results, we should stress a technical 
point. In the particular case of caustics, calculating the two func- 
tions in the same polar angle becomes a non-trivial issue. This is 
because in general, the source plane points (yit,V2t) are not equally 
distributed in angle, as they are obtained scanning angular values 
in the lens plane which map nonlinearly to angular values in the 
source plane. Thus in general, a source plane angle does not cor- 
respond to the same lens plane angle. Yet, to compute T> 2 for the 
caustics, it is necessary that both xes and xpa be calculated at the 
same polar angle position in the source plane. Thus, to enforce this 
last point, we first determine the polar angle corresponding to each 
point (yit, J/2t) obtained with the exact solution, i.e., 



fas = arctan 




and obtain the corresponding radial coordinate ies = Vt(4>s) = 
\J y\ t + y\ t . In the same way, we compute the polar angle of the 
tangential caustic obtained with the Perturbative Approach (which 
we denote by </>g A ), i.e., 

i pa , ( yit(4>h)\ 

fas = arctan x — -i , 

\yit((t)L)J 

where y lt and y 2t are given in eq. (|25} and y t PA = y/y\ t + y\ t . 
We then vary the angle far, (only the interval ^ far. ^ 7r/2 is 
needed, for symmetry reasons) such that for each radial position 
yt, the angles fas and fa% A are chosen to have the same value at 
step i. Finally, having determined (yt, fas) for the exact solution 
and (yf A , <ps A ) we proceed to compute T> 2 as in eq. J4U . 
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Fig. [2] shows T> 2 as a function of e for some values o{^ K s . 
In the left panel, the results for critical curves are shown. Since 
the perturbation increases with e, T> 2 also increases with e, as we 
might expect. In addition, T> 2 decreases as k s increases, at least for 
k s < 1.0. In the right panel, we show the results for caustics. The 
behaviour of T> 2 is qualitatively similar to that of critical curves, 
except for at the highest k s , where the behaviours are reversed. 
However, the values of T> 2 computed for caustics are higher than 
the corresponding ones for critical curves, for a given (n s ,e). This 
means that imposing cut-off values of T> 2 for matching caustics, 
we will match the corresponding critical curves automatically. We 
found by visual inspection that for D 2 ^ 5 x 1CP 4 there is a very 
good match for the caustic curves. In Fig. Q] we show the values of 
T> 2 calculated for each example, demonstrating visually the validity 
of this diagnostic measure. In particular, we have checked that cut- 
off values of V 2 higher but close to our chosen limit of 5 x 1CP 4 
are not suited for matching caustic curves well. 

T o estimate the validity of the Perturbative Approach, lAlarj 
d2007t) introduced the parameter D = 3|C 3 1 (5 x) 2 , where (Sx) 2 
corresponds to the difference between the arc contours obtained 
in the perturbative approach and the Einstein radius, and C3 is the 
third-order term in the Taylor expansion of the gravitational po- 
tential (see eq. |17| >. In order for the Perturbative Approach to be 
accurate, D should be small. For models based on the SIS profile, 
this condition is always true, since C3 = (which is consistent 
with the fact that the Perturbative Method is exact in this case). For 
other pseudo-elliptical models, usually C3 7^ 0. 

Here we adapt the definition of D to be used for critical curves, 
such that Sx is now the radial deviation of these curves with respect 
to xe- We associate a unique value of D to the tangential critical 
curve, using its maximum value over this curve, which corresponds 
to 



0.8 



= 3|C 3 |max{(a; t (0) - x E ) 2 }, 



(42) 



where x± is given in eq. d24b and ^ <f> ^ 2-k. Following Alard's 
criterion (i.e. D ma x >C 1), it would be expected that the critical 
curves and caustics obtained with the Perturbative Approach would 
be close to the ones obtained in the exact case when both e and 
k s are small. We compute D ma x for the curves shown in Fig. [T] 
obtaining D max = 0.006, 0.38 and 0.55 from left to right pan- 
els. Contrary to expectations, when D ma x increases, the curves ob- 
tained with the Perturbative Approach become more similar to the 
exact solutions. Therefore, the criterion D max <§C 1 does not reflect 
the validity of the Perturbative Approach for these cases. More- 
over, -D max is not scale-invariant (i.e. D max oc r 2 , where r B is the 
length scale of the PNFW model). These considerations show that 
this measure is not well-suited to assess the limit of validity of the 
method for caustics and critical curves. This result emphasizes the 
relevance of our definition of T> 2 as a measure for the validity of 
the Perturbative Approach for critical curves and caustics. 

For the application of our criterion, we define e^x, f° r a given 
fig, as the ellipticity threshold giving T> 2 — 5 x 10 4 . This will be 
used as a measure of the limit of applicability for the Perturbative 
Approach for critical curves and caustics. Fig. [3] shows the maxi- 
mum values of e as a function of k s for the PNFW model, for some 
cut-off values of T> 2 . The emtx( K s) function shown in this figure is 
well-fitted by a Pade approximant of the form 



4 Throughout this work, following [Dumet-Montoya e t alj d2012l) , we will 
consider the range k s ^ 1.5. 



£0.4 



T 



© 2 = 1 x 10" 3 
© 2 = 7.5 x 10" 4 
© 2 = 5 x 10" 4 




Figure 3. Maximum values of e obtained from some cut-off values of T> 2 
on caustics for the PNFW model. 



PA _ J2 n =0 a n( K s) 
S-max ^ — ,2 



(43) 



with a n = {-0.018, 0.235, -0.415, 0.565, -0.264} and b n 
{2.243, -3.709, 1.725}. 



4.2 Comparison between Deformation Cross Sections 

In this section, we compare the exact and perturbative solutions for 
the deformation cross section in order to establish limits of validity 
for the approximation of this quantity. We then contrast these limits 
to those obtained for caustics and critical curves as done in Sec. 14. 11 
(i.e. by imposing e < Emtx f° r eacn K s). If within this regime the 
Perturbative Approach and the exact solution of the deformation 
cross section do not agree well, this can impose additional limits to 
the applicability of the Perturbative Approach. 

To quantify the difference between the deformation cross sec- 
tions, we compute their relative difference 



OES,R th - TPA.Rth 



°ES,R th 



(44) 



where the subscripts ES and PA refer to the exact and perturbative 
calculations, respectively. 

In Fig.|4]we show Ao-io/Sio as a function of k s for some val- 
ues of e. In the left panel we compare the exact solution with the ex- 
pansion for low ellipticities in the Perturbative Approach, eq. 1371 , 
while in the right panel we compare with the general expression, 
eq. l |3lFl The perturbative calculation for the axially symmetric 
NFW model (e = 0, eq. |32t is a good approximation in this case, 
since Acrio/ffio < 10% for the entire allowed range of k s . For 



5 It should be noted that, due to the quadratic scaling with respect to e 
of the deviation from the axially symmetric case, eq. {37} is an excellent 
approximation for eq. Bit for low ellipticities. Thus the left panel of Fig. [4] 
would remain essentially unchanged if we used eq. )31| instead of eq. {37} 
there. 
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Figure 4. Relative deviation among deformation cross sections for the PNFW model as a function of re s for some values of e. Left panel: Actio / &10 between 
exact solution and eq. )37t . Right panel: Actio /<rio between exact solution and eq. )3U . 



values of e < 0.1 the Perturbative Approach is a good approxima- 
tion only for n s > 0.5. As e increases, the difference is larger at 
smaller values of k s . However, the perturbative calculation is accu- 
rate to within about 10% for k 3 > 0.7 up to e = 0.3 (see the right 
panel of Fig. [4]). 

Additionally, we computed Aa" t h/5"th as a function of the 
threshold Rth. We find that Aath/5'th can exceed 50% at values 
of i? t h < 2.5, since for these values of R t h, the constant distortion 
curves are far from the tangential critical curves, meaning that the 
premises of the Perturbative Approach do not apply. However, as 
i?th increases, the relative deviations among the deformation cross 
sections decrease. In particular, we found that for n s > 0.9 and 
-Rth > 7.5, these relative deviations do not depend on -Ra- 
in Fig. [5] we show isocontours of Actio /5"io, for the ex- 
act and perturbative calculations, in the parameter space n s -e to- 
gether with the curve £^ x (k s )- We see that the constraints im- 
posed by Actio /o"io and e^tx are complementary, meaning that for 
k s < 1.0 the constraint obtained with caustics and critical curves is 
the strongest, while the opposite is true for k 3 > 1.0 if we impose 
that the maximum fractional deviation for the cross section is 10%. 

We may then combine the constraints to define a region lim- 
ited approximately by the curves 

£ = | ^tx(« S ), <: 1.0, 
[ 0.33, k s > 1.0. 

Within this region the Perturbative Approach can replace the ex- 
act computation of critical curves, caustics, and deformation cross 
section with high accuracy. In other words, in this region, the Per- 
turbative Approach can reproduce well the exact local mapping of 
the lens equation. 



5 CONCLUDING REMARKS 

The Perturbative Approach dAlardll2007l I2008L 120091 boidl) pro- 
vides analytical solutions for gravitational arcs by solving the lens 
equation linearised around the Einstein Ring solution. This method 
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Figure 5. Comparison between deformation cross sections for the PNFW 
model space parameter. Contours of constant Actio /ctio. The solid line 
shows the £m^ x («: s ) curve. 



has a wide range of potential applications, from the inverse prob- 
lem in strong lensing to fast arc simulations. This technique goes 
beyond other analytical approximations in the literature in that it 
may be used for generic lens models (including mass distributions 
arising from N-body simulations) and for finite sources. 

A key aspect for practical applications of the method that has 
not been systematically addressed before is the determination of its 
limit of validity. Motivated by this issue, in this paper we aimed 
to determine the accuracy of the Perturbative Approach for caus- 
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tics and critical curves and for the deformation arc cross section. 
Although these quantities do not involve arcs (i.e. the lensing of 
finite sources) they allow one to obtain limits on the accuracy of 
the linearized mapping from the Perturbative Approach. Also, the 
parameter space to be probed is significantly decreased, since these 
quantities depend basically on the lens properties and not the source 
ones. 

We have considered a restricted set of lens models, more 
specifically those with elliptical lens potentials, and in particular 
the PNFW and SIEP models, which are nevertheless widely used 
in strong lensing applications, specially for the inverse modelling. 
Whenever possible, we sought to derive analytic expressions for 
the quantities involved in the calculations, many of which are new. 
Some are valid for the Perturbative Approach in general, others ap- 
ply to pseudo-elliptic lens models. The main results of the paper 
are summarized below. 

We obtained analytic expressions for the constant distortion 
curves in the Perturbative Approach (eqs. [27] and 1281 . which, in 
the lens plane, are found to be self-similar to the tangential critical 
curve. We derived an analytic formula for the deformation cross 
section (eg. 13 lb. which reproduces the scaling of the arc cross sec- 
tion with iith obtained numerically in previous works. For axi- 
ally symmetric models the cross section is obtained in closed form 
(eq.ES. 

We have obtained simple analytic expressions for the per- 
turbative functions for pseudo-elliptical models, which are valid 
for any choice of the ellipticity para metriza t ion (eq. |34> . These 
expressions generalize thos e given in lAlardl d2007l 120081) and in 
lHabara & Yamamotol j2(5TTh . 

We derive approximate solutions to the tangential critical 
curve (eq. [36} and for the deformation cross section (eq. 1371 for 
low ellipticities in pseudo-elliptical models. We show that the de- 
viation of the cross section with respect to the axially symmetric 
case is quadratic in the ellipticity, as noted numerically from earlier 
work. 

We have considered the SIEP and the PNFW models to rep- 
resent lenses at galaxy and galaxy cluster mass scales. We have 
shown that the Perturbative Approach provides the exact solution 
for the SIEP model. For the PNFW model, we compared the critical 
curves and caustics obtained with this approach with those obtained 
with the exact solution for a wide range of values of k 3 and e. 

W e show that the criterion D max <C 1 proposed by lAlardl 
J2007I) extended to be applied the tangential critical curve (eq.!42b. 
is not adequate to set a limit of validity for these cases. To this 
end, we use a figure-of-merit, T> 2 (eq.l4U to quantify the deviation 
of the Perturbative Approach from the exact solution for caustics 
and critical curves. We verify that T> 2 provides a quantitative de- 
scription of the deviation among both solutions. In particular, T> 2 
decreases with k 3 (as can be drawn from Fig.[T) and increases with 
e (as expected from the increasing of the perturbation to the lensing 
potential with e). Since the deviation between the exact and pertur- 
bative solutions for caustics is higher than the deviation for critical 
curves, it is sufficient to set a limit on T> 2 for caustics to ensure a 
small deviation for critical curves. 

By setting a threshold on T> 2 computed at caustics, a maxi- 
mum value of e is determined for each k s , such that a good match- 
ing for caustics and also for critical curves is ensured. We determine 
these maximum values £max( K s) by choosing V 2 = 5xl0~ 4 .This 
defines a domain of applicability of the Perturbative Approach for 
the PNFW model in the range of k s being considered. We provide 
a fitting function for e„iax( K s) (eq.!43t. For k s < 0.8, the Pertur- 



bative Approach is limited to e < 0.1. However, for k s > 1.0 it is 
possible to use this approach even up to e — 0.4 for these cases. 

Another limit on the PFNW model parameters is obtained 
from the comparison of the deformation cross section for both ex- 
act and perturbative calculations. The fractional deviation is less 
than 10% (Fig. [5} for n s > 0.7 and e < 0.3 (corresponding to 
e E < 0.55). 

We may use these results to set further constraints on the el- 
lipticity parameter of the PNFW model, by requiring an agreement 
with the exact o"K th , besides the condition e < e|^.(k s ). This en- 
sures that caustics, critical curves, and the local mapping are well 
reproduced by the Perturbative Approach for the PNFW model. 
The combined restriction, imposing the matching for caustics and 
an agreement to about 10% for deformation cross sections, is given 
in eq. <45t . 

In this paper we provided a first systematic attempt to set lim- 
its on the domain of applicability of the Perturbative Approach for 
strong lensing in terms of the parameters of a given lens model, 
more specifically for the PNFW model. The limits are imposed 
so that the caustics, critical curves and deformation cross section 
match the exact solutions with a given accuracy. Although these 
quantities are useful for strong lensing applications, it is important 
to d etermine a domain of validity for arcs/finite sources. For exam- 
ple, [Habara^&^mamot^ d201 ll) investigated the domain of valid- 
ity of the Perturbative Approach for extended circular sources. It 
is argued that Perturbative Approach can be used for sources with 
radius < 0.2 ie up to e ~ 0.3. This result should be extended for 
generic configurations probing the space of the source and lens pa- 
rameters and their relative position. We expect that the limits here 
obtained can be connected to the domain of validity for arcs provid- 
ing guidance to the exploration of this wider parameter space. The 
systematic application to arcs and connection to the current results 
is left for a subsequent work. It is also important to check whether 
the criterion established here for the T> 2 threshold can be applied 
to other lens models, so that we have an a priori criterium for the 
domain of validity of the Perturbative Approach regardless of the 
specific model. 

The usefulness of the Perturbative Approach justifies the 
search for a determination its accuracy and limit of applicability. 
Once this is established we will be able to safely use this promising 
technique in a number of applications, within its domain of validity. 
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